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The nuclear equation of state (EoS) is investigated by flow phenomena 
in relativistic heavy-ion collisions, both in transverse and radial direction, 
in comparison to experimental data from 150 A MeV to 11 A GeV. To this 
aim the collective dynamics of the nucleus-nucleus collision is described 
within a transport model of the coupled channel RBUU type. There are 
two factors which dominantly determine the baryon flow at these ener- 
gies: the momentum dependence of the scalar [Us) and vector potentials 
(C/^t) for baryons and the resonance/string degrees of freedom for ener- 
getic hadron excitations. We fix the explicit momentum dependence of the 
nucleon-meson couplings by the nucleon optical potential up to 1 GeV and 
extrapolate to higher energy. When assuming the optical potential to van- 
ish identically for Ekin > 3.5 GeV we simultaneously describe the sideward 
flow data of the PLASTIC BALL, FOPI, EoS and E877 collaborations, 
the elliptic flow data of the E895 and E877 collaborations and approxi- 
mately the rapidity and transverse mass distribution of protons at AGS 
energies without employing any explicit assumption on a phase transition 
in the EoS. However, the gradual change from hadronic to string degrees of 
freedom with increasing bombarding energy can be viewed as a transition 
from hadronic to string matter, i.e. a dissolution of hadrons at high energy 
density. 

PACS numbers: 21.65. -hf; 25.75.-q; 25.75.Ld 



1. Introduction 

Relativistic heavy-ion collisions (RHIC) provide a unique tool to study 
nuclear matter at high densities and temperatures, reminiscent of the early 
big-bang of the universe, but with better statistics and under controlled 
conditions. These reactions also provide constraints on the interior of neu- 
tron stars, where the nuclear equation-of-state (EoS) plays an essential role 
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for the possible existence of an inner quark core or an extended mixed phase 
of quarks and hadrons |jl|. However, since in a RHIC the system initially is 
far away from thermal and chemical equilibrium, both particle production 
and collective motion depend on various quantities such as the stiffness of 
the EoS, the momentum dependence of the interaction or mean- field po- 
tentials (MDI), in-medium modifications of the A^A^ cross section ctnn, the 
initial momentum distribution of the nucleons f^, ^ as well as the number 
of hadronic degress of freedom accounted for in the transport simulation 
[^. It is thus necessary not to focus on a single observable alone but to 
investigate the dynamical evolution of the RHIC within a single model that 
is able to describe all relevant single-particle as well as collective quantities. 

Whereas the experimental and theoretical studies of collective nuclear 
flow have been restricted to the 1-2 A GeV energy regime in the past 
more recently both the directed transverse flow (sideward flow) and the flow 
tensor (elliptic flow) have been measured and reported by the BNL-E877 
collaboration ||^, |^, ^, [l^ for heavy-ion {Au + Au) collisions at AGS energies 
in the energy range of 1 A GeV < Einc < 11^ GeV. In this energy range 
the directed transverse flow first grows, saturates at around 2 A GeV, and 
then decreases experimentally with energy showing no minimum as expected 
from hydro dynamical calculations including a first order phase transition in 
the EoS [|ll| . Whether this decrease in directed flow or the change of sign 
in elliptic flow is indicative of a phase transition [^] is a question of high 
current interest. 

In this contribution the collective behaviour of nuclear matter in a heavy- 
ion collision is reviewed in the energy range from 150 A MeV to 11 A GeV 
for various systems using the transport model [|l3|. For energies above 1 A 
GeV it has been, furthermore, complemented by the string dynamics from 
the HSD transport approach [14| which has been tested extensively for p+A 
and A + A collisions from SIS to SPS energies [15|. 



2. The extended RBUU-Model 



To describe the heavy-ion collision data at energies starting from the SIS 
at GSI to the SPS regime at CERN, relativistic transport models have been 
extensively used |14, |l8[ 19 1. For a general derivation of transport 

theories the reader is refered to Ref. [20| and to Ref. for a recent 

review. Among these transport models the Relativistic Boltzmann-Uehling- 
Uhlenbeck (RBUU) approach incorporates the relativistic mean-field (RMF) 
theory, which is applicable also to various nuclear structure problems as well 
as for neutron star studies [Q, Since it is based on an effective hadronic 
Lagrangian density it allows to evaluate directly the nuclear EoS at zero and 
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finite nuclear temperature T as well as the scalar and vector mean fields Us 
and C/^, that determine the in- medium particle properties, for arbitrary 
configurations of nucleons in phase space [EOl . Here we essentially base the 



studies on the Lagrangian (parameter set NL3) from Ref. |22| since this 
Lagrangian has been applied widely in the analysis of heavy-ion collisions 



by various groups [O, 23, p4 



We recall that the most simple versions of RMF theories assume the 
scalar and vector fields to be represented by point-like meson-baryon cou- 
plings. These couplings lead to a linearly growing Schrodinger-equivalent 
potential in nuclear matter as a function of the kinetic energy -E^m! which 
naturally explains the energy dependence of the nucleon optical potential 
at low energies (< 200 MeV). However, a simple RMF does not describe 
the nucleon optical potential at higher energies, where it deviates substan- 
tially from a linear function and saturates at -Efci^ ~ 1 GeV |2^. Since the 
energy dependence of sideward flow is controlled in part by the nucleon op- 
tical potential, the simple RMF cannot be applied to high-energy heavy-ion 
collisions. In order to remedy this aspect, the more sophisticated RBUU ap- 
proaches invoke an explicit momentum dependence of the coupling constant. 



i.e. a form factor for the meson-baryon couplings |14, [T^, 26|. 

Further important ingredients at AGS energies are the resonance/string 
degrees of freedom which are excited during the reaction in high energy 
baryon-baryon or meson-baryon collisions. While at SIS energies particle 
production mainly occurs through baryon resonance production and their 
decay, the string phenomenology is found to work well at SPS energies (~ 
200 A GeV) jl^. One of the characteristic features of the AGS energy 
regime is the competition between these two particle production mecha- 
nisms which might be separated by some energy scale y/Sgw 1 27]. Due to 
this complexity there are various ways to implement elementary cross sec- 



tions in transport models M, 12, 15, u% ESf. One of the extremes is 



to parameterize all possible cross sections directly for multi-pion produc- 
tion, A'^A^ — > NNnir (n > 3) only through N, A, vr degrees of freedom; the 
other extreme is to fully apply string phenomenology in this energy region 
without employing any resonances. Although it is possible to reproduce 
the elementary cross sections from NN and ttN collisions and the inclusive 
final hadron spectra in heavy-ion collisions within these different models, 
we expect that differences should appear in the dynamical evolution of the 
system, e.g. in the thermodynamical properties Q ^] and in collective 
flow |2^. For example, if thermal equilibrium is achieved at a given energy 
density, models with a larger number of degrees of freedom including strings 
will give smaller temperature and pressure 

We have employed here the combination of a resonance production model 
[T| and the Lund string model |32] as incorporated in the Hadron-String- 



4 



TEXT PRINTED ON FEBRUARY 8, 2008 



Dynamics (HSD) approach |14, 15|. In the practical implementation for 
NN or MN collisions at invariant energies lower (higher) than a threshold 
energy ^/Sgw resonances (strings) are assumed to be excited (see below). 



2.1. The optical potential 



In this presentation the scalar and vector mean fields Us and C/^ are 
calculated on the basis of the same Lagrangian density as considered in 
Ref. 1 13], which contains nucleon, a and lo meson fields and nonlinear self- 
interactions of the scalar field (cf. NL3 parameter set 123] ). The scalar and 



vector form factors at the vertices are taken into account in the form [14| 



fsip) 



A 



A2 +p2 



and 



Up) 



A2 +p2 



(1) 



where the cut-off parameters A^ = 1.0 GeV and A^ = 0.9 GeV are obtained 
by fitting the Schrodinger equivalent potential, 



U sep{,Ekin) 



to Dirac phenomenology for intermediate energy proton-nucleus scattering 




EKin(GeV) 



Fig. 1 The Schrodinger equivalent potential (g) at density po as a function of the 
nucleon kinetic energy Ekin ■ The solid curve (RBUU) results from the momentum- 
dependent potentials discussed in the text. The data points are from Hama et al. 

Ill- 

The resulting Schrodinger equivalent potential (^) is shown in Fig. 1 as 
a function of the nucleon kinetic energy with respect to the nuclear matter 
at rest in comparison to the data from Hama et al. ||2^ (open squares). The 
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experimental increase of the Schrodinger equivalent potential up to E^in = 1 
GeV is decribed quite well; then the potential decreases and is set to zero 
above 3.5 GeV. 

For the transition rate in the collision term of the transport model we 
employ in-medium cross sections as in Ref. that are parameterized in 
line with the corresponding experimental data for ^/s < yjsgw For higher 
invariant collision energies we adopt the Lund string formation and frag- 
mentation model as incorporated in the HSD transport approach |14| 
which has been used extensively for the description of particle production in 
nucleus-nucleus collisions from SIS to SPS energies ||l^ . In the present rela- 
tivistic transport approach (RBUU) as in Ref. |13] we explicitely propagate 
nucleons and A's as well as all baryon resonances up to a mass of 2 GeV 

1. Furthermore, 



with their isospin degrees of freedom |31 



vr,??, p, to, 



K,K 



and a mesons are propagated, too, where the a is a short lived effective 
resonance that describes s-wave tttt scattering. For more details we refer 
the reader to Refs. 33 1 concerning the low energy cross sections and to 
Refs. 1 14, 15 1 with respect to the implementation of the string dynamics. 



2.2. Transverse mass spectra of protons 

In Fig. 2 we show the dependence of the calculated proton transverse 
mass spectra in a central collision of Au + Au at 11.6 A GeV for b < 3.5 
fm for ^/Ssw = 2.6 GeV (dotted histograms) and 3.5 GeV (solid histograms) 
in comparison to the experimental data of the E802 collaboration ||3^. A 
cascade calculation (crosses) is shown additionally for ^/sJ^ = 3.5 GeV to 
demonstrate the effect of the mean-field potentials which lead to a reduction 
of the transverse mass spectra below 0.3 GeV. As expected, the transverse 
mass spectrum is softer for smaller ^/s^ due to the larger number of degrees 
of freedom in the string model relative to the resonance model. 
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Fig. 2: The transverse mass spectra of protons for Au + Au collisions at 6 < 3.5 
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fm. The solid line and the dot- dashed line with crosses are results for ^Ssw—^-^ 
GeV with and without nuclear potentials, respectively. The dotted line RBUU(2.6) 
is for GeV. The data points are taken from the E802 collaboration js^. 

We note that strings may be regarded as hadronic excitations in the 
continuum of lifetime tp ~ 0.8 fm/c (in their rest frame) that take over 
a significant part of the incident collision energy by their invariant mass. 
They decay dominantly to light baryons and mesons and only to a low 
extent to heavy baryon resonances. Thus the number of particles for fixed 
system time is larger for string excitations than for the resonance model 
where several hadrons propagate as a single heavy resonance which might 
be regarded as a cluster of a nucleon + n pions. As a consequence the 
translational energies are suppressed in string excitations and, as a result, 
the temperature as well as the pressure are smaller when exciting strings. 

From the above comparison with the experimental transverse mass spec- 
tra for protons we find ^/s^ ~ 3.5 GeV, which implies that binary final 
baryon channels should dominate up to ~ 3.5 GeV which corresponds 
to a proton laboratory energy of about 4.6 GeV for pp collisions. 



3. SIS energies 

3.1. Transverse flow 

Within the RBUU model described above we now calculate the trans- 
verse flow for various systems and beam energies and analyse the dependence 
on different quantities. The flow F is defined as the slope of the transverse 
momentum distribution at midrapidity, 

F='-^ , (3) 

which is essentially generated by the participating matter in the 'fireball' Q. 
The latter finding explains why flow does not clearly distinguish between 
an EoS with and without momentum-dependent forces. Since the fireball 
contains the stopped matter, the relative momenta in the fireball (besides 
the unordered thermal motion) are small. Only when applying additional 
cuts, e.g. on high transverse momenta i.e. by selecting particles 

escaping early from the fireball, or selecting mainly participant or spectator 
particles by appropriate 0cm-Cuts [^5| , a difference between the momentum- 
dependent and momentum-independent EoS can be established. 
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The FOPI data on proton flow Q indicate a decrease of sidewards flow 
above 1 A GeV incident energy following the well known logarithmic in- 
crease at low energies. Using standard potential parameterizations, both 
nonrelativistic Q and relativistic |13], this behavior cannot be understood 
within conventional transport models. In the latter the optical potential 
stays constant or even increases at high momenta and therefore the repul- 
sion generated from the momentum-dependent forces in a HIC gives rise to 
a significant contribution to the flow signal. However, since the nucleon- 



nucleus optical potential is only known up to 1 GeV experimentally |25|, it 
was proposed that this decrease of flow above 1 A GeV might indicate a 
decrease of the optical potential at high relative momenta or at high baryon 
density p^ . 

In Ref. Q the transverse momentum of the baryons has been disentan- 
gled into a collisional part, a mean- field part and a part originating from 
the Fermi-motion of the particles. 



Pt 



Pt 



Fermi 



(4) 



where p^*^''™' practically does not contribute at midrapidity. Thus disentan- 
gling the flow signal into a collisional and a potential part, it turns out that 
~50 % of the flow stems from the particle - particle collisions while roughly 
another 50 % are generated by the potential repulsion at 1 A GeV. Fig. 
3 shows these two contributions for flow for the system Ni + Ni at b=4 
fm in comparison to the experimental data using different EoS denoted by 
'hard', 'medium' and 'soft'. Both the collisional "background" and the po- 
tential part rise up to 1 A GeV incident energy and remain constant above, 
whereas the data indicate a decrease above 1 A GeV. As seen from Fig. 3 
the Ni + Ni data are described reasonably by both a 'soft' and 'medium' 
EoS, while a cascade calculation fails substantially. 
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Fig. 3: Transverse flow for a Ni + Ni collision at b=4 fm as calculated within 
the RBUU model in cascade mode (dashed lower line) and for equations of state 
with different compressibilities K in comparison to data from EOS, Plastic Ball 
and FOPI as compiled by Ref. . 

A note of caution has to be added here: The flow F @ not only depends 
on the baryon self-energies Us and but also on the number of (resonance) 
degrees of freedom above about 1 A GeV as first pointed out by Hombach 
et al. This observation will become crucial at AGS energies. 



3.2. Radial flow 



Radial flow has been discovered [36| when analyzing the flow pattern of 
very central events of RHIC. In contrast to transverse flow up to about 70 % 
of the incident energy (stored in the hot compressed fireball) is released as 
ordered radial expansion of the nuclear matter. Thus the hope is to extract 
information especially on the compressibility of the EoS via the magnitude 
of the radial flow. 

Experimentally the radial flow is characterized or fitted in terms of the 
Siemens-Rasmussen formula [BtII 



p ■ e 



-7E/T 



f sinh a 



V 



a 



{jE + T) - T ■ cosha^ 



(5) 



with 7 = (1 — /5^)~"'^/^ and a = j(3p/T, while /? denotes the flow velocity 
and T characterizes some temperature. We follow the same strategy in 
our RBUU calculations and apply a least square fit to the RBUU nucleon 
spectra using Eq. dsl). 
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Fig. 4: The temperature (l.h.s.) and radial flow velocity (r.h.s.) for central 
Au + Au collisions evaluated via Eq. (|) from the RBUU calculations in comparison 
to the experimental data from Refs. |3q, P9^, M. The symbol 's' denotes a soft EoS 
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without momentum dependent forces, 'h' a hard EoS and 'smd', 'mmd' and 'hmd' 
correspond to a soft, medium and hard momentum dependent EoS, respectively. 



The results of the RBUU calculations for central Au + Au collisions are 
shown in Fig. 4 as a function of the bombarding energy in comparison to 



the data from [ |38| , |39| , 4C|. We find that the 'temperature' T is systemat- 
ically underpredicted in all schemes investigated (soft and hard EoS, with 
and without momentum dependent forces), and that the flow velocities are 
not correctly reproduced, being too low at low energy and crossing the ex- 
perimental data around 800 A MeV. This might indicate a strong binding 
from the potential which gives not enough repulsion at high densities and 
overcompensates the collisional pressure from the fireball. However, the 
nucleon spectra resulting from the RBUU calculation show a strong non- 
thermal component at low incident energies and are thus in contradiction to 
the physical picture behind Eq. (jsl) which assumes an isentropic expansion 



of a thermal equilibrated source. In Ref. |41| the degree of equilibration in 
a HIC has been investigated as a function of the incident energy and the 
system mass and it has been found that even the most massive systems like 
Au + Au do not equilibrate at low energies. 



4. AGS energies 

4..I. Sidewards flow 

We now turn to the AGS energy regime from 1-11 AGeV. The calcula- 
tions are performed for the impact parameter h = 6fm for Au + Au systems, 
since for this impact parameter we get the maximum flow which corresponds 
to the multiplicity bins M3 and M4 as deflned by the Plastic Ball collabo- 
ration [^] at BEVALAC/SIS energies. In Fig. 5 (l.h.s.) the transverse flow 
(^) is displayed in comparison to the data from Refs. ^, for Au + Au 
systems. The solid line (RBUU with y/Ssw = 3.5 GeV) is obtained with 
the scalar and vector self energies as discussed above, Eq. (|l|). The dotted 
line (CASCADE with ^Ssw = 3.5 GeV) corresponds to cascade calculations 
for reference in order to show the effect of the mean field relative to that 
from collisions. We observe that the solid line (RBUU, cf. Fig. 1) is in good 
agreement with the flow data at all energies; above bombarding energies 
of 6 A GeV the results are practically indentical to the cascade calculations 
showing the potential effects to cancel out. 

We note that the sideward flow shows a maximum around 2 A GeV for 
Au + Au and decreases continuously at higher beam energy (> 2 A GeV) 



10 



TEXT PRINTED ON FEBRUARY 8, 2008 



without showing any exphcit minimum as in Ref. |11]. This is due to the 
fact that the repulsive force caused by the vector mean field decreases at 
high beam energies (cf. Fig. 1) such that in the initial phase of the collision 
there are no longer strong gradients of the potential within the reaction 
plane. In subsequent collisions, which are important for Au + Au due to 
the system size, the kinetic energy of the particles relative to the local rest 
frame is then in a range (-Bfcm < 1 GeV) where the Schrodinger equivalent 
potential (at density po) is determined by the experimental data [^]. We 
thus conclude that for the sideward flow data up to 11 A GeV one needs 
a considerably strong vector potential at low energy and that one has to 
reduce the vector mean field at high beam energy in line with Fig. 1. In 
other words, there is only a weak repulsive force at high relative momenta 
or high densities. 
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Fig. 5: (l.h.s.) The sideward flow F{y) as a function of the beam energy per 
nucleon for Au + Au collisions at & = 6 fm from the RBUU calculations. The solid 
line results for the parameter set RBUU, the dotted line for a cascade calculation 
with y^Ssw = 3.5 GeV. The data points are from the FOPI and EoS Collaborations 
[0i ill- (r.h.s.) The elliptic flow V2 of protons versus the beam energy per nucleon 
for Au + Au collisions at 5 = 6 fm from the RBUU calculations. The solid line 
results for the parameter set RBUU, the dotted line for a cascade calculation with 
y^Ssw = 3.5 GeV. The data points are from the EoS Collaboration [lo| . 

Another aspect of the decreasing sideward flow can be related to the 
dynamical change in the resonance/string degrees of freedom as already 
discussed above. For instance, for ^/Ssw = 2.6 GeV the calculated flow 
turns out to be smaller than the data above 1.5 A GeV and approaches 
the cascade limit already for ~ 3-4 A GeV. This is due to the fact that in 
strings the incident energy is stored to a larger extent in their masses and 
the translational energy is reduced accordingly. 



4-2. Elliptic flow 
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Apart from the in-plane flow of protons the out-of-plane collective flow 
provides additional information and constraints on the nuclear potentials 
involved. In this respect the elliptic flow for protons 

V2 = {{Pl-P^)/{Pl + P^)) (6) 

for \y/ycm\ ^ 0.2 is shown in Fig. 5 (r.h.s.) as a function of incident energy 
for Au + Au collisions at 6 = 6 fm. The solid line (RBUU with yjssw = 
3.5 GeV) is obtained with the same mean fields as discussed before while 
the dotted line (CASCADE with yjsgw = 3.5 GeV) stands again for the 
cascade results. The flow parameter V2 changes its sign from negative at 
low energies (< 5 A GeV) to positive elliptic flow at high energies (> bA 
GeV). 

This can be understood as follows: At low energies the squeeze-out of 
nuclear matter leads to a negative elliptic flow since projectile and target 
spectators distort the collective expansion of the 'fireball' in the reaction 
plane. At high energies the projectile and target spectators do not hinder 
anymore the in-plane expansion of the 'fireball' due to their high velocity 

c); the elliptic flow then is positive. The competition between squeeze- 
out and in-plane elliptic flow at AGS energies depends on the nature of the 



nuclear force as pointed out already by Danielewicz et al. |12|. We note, 
however, that in our calculation with the momentum-dependent potential 
(Fig. 1) we can describe both the sideward as well as elliptic flow data [^, |l^ 
simultaneously without incorporating any phase transition in the EOS as in 



Ref. PI 



In the cascade calculation the elliptic flow from squeeze-out is weaker 
due to the lack of a nuclear force which demonstrates the relative role of 
the momentum-dependent nuclear forces on the V2 observable below bom- 
barding energies of about 5 A GeV. 



5. Summary 

In this contribution we have explored the dependence of transverse and 
radial flow signals on various model inputs - that are related to the nuclear 
EoS - using the coupled channel RBUU model. We find that the mass 
distribution of the resonances included in the model plays an important 
role for the description of transverse fiow above 1 A GeV. For the radial 
flow we have concentrated on the difference between the results for the flow 
temperature T and flow velocity (3 when using different EoS. However, no 
sizeable sensitivity to the compressibility of the EoS could be established. 

On the other hand, we found that in order to reach a consistent under- 
standing of the nucleon optical potential up to 1 GeV, the transverse mass 
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distributions of protons at AGS energies as well as the excitation function 
of sidewards and elliptic flow up to 11 A GeV, the strength of the 

vector potential has to be reduced in the RBUU model at high relative 
momenta and/or densities. Otherwise, too much flow is generated in the 
early stages of the reaction and cannot be reduced at later phases where the 
Schrodinger equivalent potential is experimentally known. This constrains 
the parameterizations of the explicit momentum dependence of the vector 
and scalar mean fields and Us at high relative momenta. 

In addition, we have shown the relative role of resonance and string de- 
grees of freedom at AGS energies. By reducing the number of degrees of 
freedom via high mass resonances one can build up a higher pressure and/or 
temperature of the 'fireball' which shows up in the transverse mass spectra 
of protons as well as in the sidewards fiow [27|. A possible transition from 
resonance to string degrees of freedom is indicated by the RBUU calcula- 
tions at invariant baryon-baryon collision energies of ^/s ~ 3.5 GeV which 
corresponds to a proton laboratory energy of about 4.6 GeV. Due to Fermi 
motion of the nucleons in Au + Au collisions the transition from resonance 
to string degrees of freedom becomes smooth and starts from about 3 A 
GeV; at 11 A GeV practically all initial baryon-baryon collisions end up in 
strings, i.e. hadronic excitations in the continuum that decay to hadrons 
on a time scale of about 0.8 fm/c in their rest frame. This initial high- 
density string matter (up to 10 po at 11 j4 GeV) should not be interpreted 
as hadronic matter since it implies roughly 5 constituent quarks per fm^, 
which is more than the average quark density in a nucleon. 

It is interesting to note that at roughly 4 A GeV in central Au + Au 
collisions the ratio /t:'^ is enhanced experimentally relative to transport 
calculations [Q] . One might speculate that a restoration of chiral symmetry 
could be responsible for the softer collective response as well as enhanced 
strangeness fraction. 
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